AD-A207  829 


*nc  fiif  mpv  , 


Naval  Environmental  Prediction  Research  Facility 

Monterey,  CA  93943-5006 

Contractor  Report  CR  89-03  May  1989 


© 


SATELLITE-DERIVED 
MOISTURE-BOGUSING  PROFILES 

FOR  THE 

NORTH  ATLANTIC  OCEAN 


Lee  W.  Eddington 

Range  Operations  Department 
Pacific  Missile  Test  Center 
Point  Mugu,  CA  93042-5000 


DT1C 

ELECTE I 
MAY  0  9 1989 


approved  for  public  release,  distribution  is  unlimited 

89  5  09  059 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  This  PAGE 


la  REPORT  SECURITY  CLASSIFICATION 

UNCLASSIFIED 


2a.  SECURITY  CLASSIFICATION  AUTHORITY 


REPORT  DOCUMENTATION  PAGE 


'b  RESTRICTIVE  MARKINGS 


2b.  DECLASSIFICATION  /  DOWNGRADING  SCHEDULE 


4.  PERFORMING  ORGANIZATION  REPORT  NUMBER(S) 

GSTN-145 


2  DISTRIBUTION/ AVAILABILITY  OF  REPORT 

Approved  for  public  release; 
distribution  is  unlimited 


S  MONITORING  ORGANIZATION  REPORT  NUMBER(S) 

CR  89-03 


6a  NAME  OF  PERFORMING  ORGANIZATION  6b  OFFICE  SYMBOL  7a.  NAME  OF  MONITORING  ORGANIZATION 

Geophys.  Div. ,  Range  Ops.  Dept. ,  (if  applicable)  Naval  Environmental  Prediction 

Pacific  Missile  Test  Center  Code  3200  Research  Facility 


be  AODRESS  (City,  Stitt,  and  ZIP  Code) 


7b  ADDRESS  (City,  State,  and  ZIP  Code) 


Point  Mugu,  CA  93042-5000  gj  'b'  q 


Monterey,  CA  93943-5006 


8a.  NAME  OF  FUNDING  /  SPONSORING 
ORGANIZATION 

Office  of  Naval  Research 


8c.  ADDRESS  (City,  State,  and  ZIP  Code) 

Arlington,  VA  22217 


8b  OFFICE  SYMBOL  9.  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 
(If  applicable) 


0NT,  Code  22 


10  SOURCE  OF  FUNDING  NUMBERS 


PROGRAM 
ELEMENT  NO 

62435N 


PROJECT 

TASK 

NO 

NO 

RM35G82 

4 

WORK  UNIT 
ACCESSION  NO 

DN658752 


1 1  TITLE  (include  Security  Classification ) 

Satellite-Derived  Moisture-Bogusing  Profiles  for  the  North  Atlantic  Ocean  (U) 


12  PERSONAL  AUTHOR(S) 


13a.  type  of  REPORT 

Final 


Eddington,  Lee  W. 


13b  TIME  COVERED  14  DATE  OF  REPORT  (Year.  Month,  Day)  IlS  PAGE  COUNT 

FROM  1/88  TO  1 2/88  1989,  May  j  34 


COSATI  COOES 


GROUP  SUB-GROUP 


02 


18.  SUBJECT  TERMS  ( Continue  on  reverse  if  necessary  and  identify  by  block  number ) 

Bogusing,  -/Satellite  meteorology,  ( .  j,  ) 

Moisture  profile^- 

Cloud  analysis  j _  _ 


ABS^R^C  y(<jo ytmue  on  reverse  if  necessary  and  identify  by  block  number) 

<V\  set  of  12  vertical  profiles  of  relative  humidity  for  use  in  satellite  bogusing  of 
moisture  into  numerical  weather  prediction  model  initial  analysis  is  presented  for  the  North 
Atlantic  Ocean  region.  The  profiles  are  based  on  a  satellite  image  cloud  classification 
scheme  applied  to  a  concurrent  data  set  of  satellite  images  and  upper-air  soundings. 

Relative  humidity  profile  composites  and  statistical  significance  tests  are  presented  for 
each  cloud  category.  Sources  of  error  in  the  compositing  process  are  discussed,  and 
possible  improvements  are  suggested.  }\€,,  i  c  • 


20  DISTRIBUTION /AVAILABILITY  OF  ABSTRACT 

B UNCLASSIFIED/UNLIMITED  □  SAME  AS  RPT  □  OTIC  USERS 


22a  NAME  OF  RESPONSIBLE  INDIVIDUAL 

Dr.  Paul  M.  Tag,  contract  monitor 


OO  FORM  1473,  84  MAR  83  APR  editi-  *i  may  b«  used  until  exhausted 

All  other  editions  are  obsolete 


21.  ABSTRACT  SECURITY  CLASSIFICATION 

UNCLASSIFIED 


22b  TELEPHONE  (>r*c/vde  Area  Code)  22c.  OFFICE  SYMBOL 

(408)  647-4737  NEPRF  WU  6.2-9 


il  exhausted  SECURITY  CLASSIFICATION  OF  This  PAGE 

i0lete  UNCLASSIFIED 


CONTENTS 


J 

1.  Introduction . 1 

t  2.  Data  and  Procedures . 2 

3.  Results  .  . . 4 

4.  Discussion  and  Suggestions  .  22 

5.  Summary . 26 

Acknowledgements  .  27 

References . 28 

Distribution  .  29 


1 


INTRODUCTION 


1  . 


Numerical  weather  prediction  models  are  highly  dependent  on 
analyses  of  initial  atmospheric  conditions  in  order  to  accurately 
predict  future  conditions.  In  general,  the  more  accurate  these 
initial  analyses,  the  more  accurate  the  prediction.  The  analyses 
are  usually  good  over  land  areas  due  to  the  relatively  dense 
networks  of  surface  and  upper-air  observation  stations.  Unfortu¬ 
nately,  over  the  oceans  there  are  very  few  of  tjiese  conventional 
observations.  Thus  the  accuracy  of  the  model  initial  analyses 
for  the  ocean  areas  is  often  questionable. 

f j 

Due  to  .this  lack  of  observational  data  over  oceans,  a  tech¬ 
nique  was  created  which  uses  meteorological  satellite  image  anal¬ 
ysis  to  better  assess  the  spatial  distribution  of  moisture  in  the 

I 

atmosphere.  The  technique,  called  'moisture  bogusing.'  has  been  t'S?  «■' 
-ut,i  1  i  zed- s  i  nee  the  early  1970's  by  the  National  Meteorological 
Center  (NMC)  to  improve  their  model  ini tial ization ,  (Chu ,  1977; 

Smigielski  et  al  .  ,  1982;  Timchalk,  1 986 )  .  ^^he  technique  is 

utri  Jri2€Tl  by  assigning  one  of^  a  number  of  previously  computed 
vertical  moisture  profiles  to  different  areas  on  a  satellite 
image,  gridding  the  analysis,  and  incorporating  these  profiles 
into  model  initialization.  The  vertical  moisture  profiles  are 
computed  by  averaging  soundings  of  relative  humidity  for  differ¬ 
ent  synoptic  scale  cloud  patterns  observed  on  satellite  imagery. 

The  importance  of  moisture  bogusing  in  the  numerical  forecasting 


« 


of  precipitation  was  demonstrated  by  Lyons  (1986a)  ,  and  modifica¬ 
tions  to  NMC ’ s  moisture  profiles  were  suggested  in  Lyons  (1986b). 

This  report  presents  the  results  of  an  ongoing  study  funded 
by  the  Naval  Enviromental  Prediction  Research  Facility  (NEPRF)  to 
establish  sets  of  moisture  bogusing  profiles  for  different  ocean 
areas.  The  results  presented  here  are  for  the  North  Atlantic 
Ocean.  The  paper  proceeds  as  follows:  Section  2  describes  the 
data  and  procedures,  Section  3  describes  the  results,  Section  4 
contains  discussion  and  suggestions,  and  Section  5  is  a  summary. 

2.  DATA  AND  PROCEDURES 

In  order  to  determine  a  set  of  vertical  moisture  profiles 
for  the  North  Atlantic,  soundings  of  relative  humidity  versus 
pressure  were  matched  with  concurrent  visual  and  infrared 
Northern  Hemisphere  satellite  image  mosaics  taken  from  the  NOAA-2 
polar  orbiting  satellite.  A  total  of  471  soundings,  taken  from 
four  ocean  weather  ships  and  one  island  station,  for  the  months 
of  January,  July,  and  October  1974  were  used  in  creating  the 
bogusing  profiles.  Table  1  gives  the  station  names,  locations, 
number  of  soundings  by  month,  and  total  number  of  soundings. 
Relative  humidity  values  were  computed  from  the  upper-air  sound¬ 
ings  of  temperature  and  dewpoint  at  mandatory  and  significant 
pressure  levels.  These  values  were  interpolated  to  19  pressure 
levels  ranging  from  1000  to  100  mb  at  50  mb  increments.  Missing 
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Table  1. 

Number  of  Observations 


Station 

i 

Position 

\ 

Jan  . 

i 

Jul  . 

t 

Oct  . 

\ 

Total 

Ship 

H 

i 

38 

N 

72 

W 

1 

40 

1 

0 

i 

0 

i 

40 

Ship 

I 

i 

60 

N 

19 

w 

! 

29 

i 

30 

! 

33 

i 

92 

Ship 

J 

1 

53 

N 

19 

w 

I 

38 

! 

36 

1 

39 

i 

113 

Ship 

K 

i 

45 

N 

16 

w 

I 

41 

i 

43 

! 

46 

! 

130 

Bermuda 

i 

32 

N 

65 

w 

l 

34 

! 

31 

I 

31 

! 

96 

! 

182 

1 

140 

! 

149 

i 

471 

i - i - i 


values  above  300  mb  were  set  to  5  percent. 

A  scheme  different  from  the  NMC  cloud  pattern  categorization 
scheme  was  created  in  order  to  simplify  the  subjective  analysis. 
This  categorization  scheme  is  presented  in  Table  2.  The  first 
11  categories  are  the  identifiable  categories  running  roughly 
from  dry  to  moist.  Category  12  is  a  composite  of  all  cases  on 
which  the  analyst  was  undecided.  Category  13  is  the  composite 
ensemble  of  all  471  soundings  used  in  this  study. 

Each  sounding  was  assigned  to  one  of  the  first  12  categories 
based  on  a  subjective  analysis  of  the  corresponding  satellite 
image(s).  Composites  and  67  percent  confidence  limits  were  com¬ 
puted  for  each  category.  Significance  tests  using  the  Wilcoxon 
rank-sum  test  were  made  to  determine  the  statistcal  significance 
of  the  individual  category's  relative  humidity  differences  from 
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Table  2. 

Cloud  Pattern  Classification  Scheme 


1 - 

!  Category  * 

!  Description 

i  1 

!  Clear 

!  2 

!  Stratus 

'  3 

!  Stratocumulus 

'  4 

!  Open-celled  cumulus 

i  5 

!  Open-celled  cumulonimbus 

'  6 

'  Altocumulus  (chaotic) 

'  7 

!  Altostratus  (organized) 

'  8 

!  Thin  cirrus 

'  9 

!  Thick  cirrus 

i  10 

t  Low,  mid,  and  high  clouds 

(chaotic)  1 

'  1  1 

!  Low,  mid,  and  high  clouds 

(organized)  ' 

'  12 

1  Undecided 

!  13 

!  Ensemble 

ensemble.  These  results  are  presented  in 

the  following 

section . 

3.  RESULTS 

The  results  of  the  compositing  are  presented  in  Table  3. 

The  composites  are  also  presented  graphically  in  Figs.  1-13  with 
67  percent  confidence  limits,  and  level  by  level  significance 
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Table  3. 

Relative  Humidity  Composites  (%) 


!  Press . 

! 

Category 

(mb) 

!  - 

— 

!  - 

— 

— 

— 

!  - 

— 

!  - 

— 

!  - 

— 

!  - 

— 

!  - 

— 

!  - 

— 

!  - 

— 

1  “ 

- I 

— 

1 

1 

! 

2 

3 

4 

I 

5 

I 

6 

I 

7 

8 

! 

9 

! 

10 

! 

1  1 

1 

12 

13 

1000 

1 

76 

! 

90 

81 

74 

! 

73 

! 

81 

I 

82 

82 

! 

84 

! 

80 

| 

85 

I 

79 

80 

950 

I 

75 

1 

81 

82 

78 

I 

75 

I 

88 

I 

90 

80 

| 

81 

I 

81 

} 

87 

* 

79 

80 

900 

1 

66 

1 

66 

72 

75 

I 

75 

! 

90 

i 

92 

72 

1 

80 

! 

79 

i 

85 

1 

73 

74 

850 

1 

49 

! 

53 

57 

65 

I 

71 

I 

71 

i 

88 

60 

| 

67 

( 

73 

| 

80 

1 

65 

64 

800 

1 

44 

! 

49 

47 

56 

f 

65 

j 

65 

I 

78 

55 

I 

61 

j 

65 

I 

79 

I 

57 

56 

750 

1 

42 

I 

42 

44 

49 

| 

61 

! 

62 

I 

69 

50 

I 

64 

! 

59 

1 

80 

I 

52 

52 

700 

1 

40 

{ 

37 

39 

41 

I 

55 

! 

63 

! 

67 

46 

! 

67 

| 

57 

I 

75 

1 

48 

48 

650 

I 

38 

f 

35 

37 

36 

1 

52 

( 

63 

I 

62 

44 

! 

70 

j 

55 

I 

75 

l 

46 

45 

600 

1 

36 

i 

31 

32 

34 

! 

51 

! 

53 

! 

64 

42 

I 

64 

! 

56 

i 

75 

i 

42 

43 

550 

I 

30 

30 

33 

32 

{ 

46 

I 

53 

I 

54 

43 

! 

70 

1 

52 

| 

\ 

71 

i 

41 

4  1 

500 

! 

24 

| 

31 

33 

29 

i 

33 

I 

46 

i 

45 

44 

i 

67 

I 

48 

l 

66 

l 

38 

39 

450 

I 

21 

I 

32 

31 

26 

! 

18 

I 

36 

! 

44 

46 

! 

62 

I 

48 

( 

I 

54 

! 

34 

36 

400 

i 

20 

1 

28 

28 

18 

1 

8 

l 

33 

I 

33 

40 

I 

49 

i 

42 

! 

42 

! 

30 

30 

350 

I 

25 

I 

26 

25 

13 

I 

5 

I 

29 

! 

24 

30 

i 

26 

I 

28 

! 

23 

1 

24 

23 

300 

i 

13 

! 

13 

9 

7 

l 

5 

i 

5 

? 

9 

9 

! 

12 

i 

8 

1 

9 

I 

8 

9 

250 

| 

5 

I 

5 

5 

5 

i 

5 

l 

5 

l 

5 

! 

5 

! 

9 

| 

5 

* 

5 

I 

5 

5 

200 

! 

5 

i 

5 

5 

5 

l 

5 

1 

5 

l 

5 

I 

5 

! 

5 

i 

5 

| 

5 

I 

5 

5 

150 

I 

5 

I 

5 

5 

5 

| 

5 

! 

5 

l 

5 

I 

5 

I 

5 

! 

5 

! 

5 

l 

5 

5 

100 

! 

5 

I 

5 

5 

5 

! 

5 

I 

5 

i 

- 1 

5 

! 

5 

! 

5 

I 

-  ! 

5 

! 

-  ! 

5 

1 

5 

5 

5 


values  which  compare  the  category  composite  to  the  rest  of  the 
sample.  Figure  1  shows  the  ensemble  (category  13)  with  67 
percent  confidence  limits  for  each  level.  The  confidence  limits 
for  this  and  all  other  categories  were  computed  using  nonpara- 
metric  statistics.  Nonparametr ic  statistics  were  used  because 
they  do  not  need  a  normally  distributed  sample  to  be  valid.  This 
choice  seemed  appropriate  when  dealing  with  relative  humidity 
values  with  limiting  values  of  0  and  100  percent  and  asymmetric 
distributions.  Confidence  limits  computed  in  this  manner  repre¬ 
sent  the  upper  and  lower  bounds  of  the  middle  two-thirds  of  the 
relative  humidity  values  observed  at  that  level. 

The  number  of  soundings  in  each  category  is  given  in  the  top 
right  hand  corner  of  the  figures  as  “N=‘  followed  by  two  numbers. 
The  first  number  is  the  total  number  of  soundings.  The  second 
number,  within  the  parenthesis,  is  the  number  of  soundings  with 
data  above  400  mb.  The  reason  for  this  difference  is  that  the 
data  base  of  soundings  from  Ship  H  and  Bermuda  had  no  data  above 
400  mb. 

Significance  values  were  only  computed  for  the  first  15 
pressure  levels  (1000-300  mb)  due  to  the  lack  of  significant 
variability  in  the  data  above  300  mb.  The  significance  values 
represent  the  probability  that  the  individual  category  versus  the 
ensemble  relative  humidity  rank-sum  discrepancies  was  due  to 
chance.  For  example,  a  value  of  0.05  (denoted  by  the  dashed  line 
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FIG.  3.  Same  as  in  Fig.  2  but  for  Category  2  (stratus) 
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PRESSURE  (MB) 


CATEGORY  3 


FIG . 


Same  as  in  Fig.  2  but  for  Category  3  (stratocumulus) . 


PRESSURE  (MB) 


FIG.  5.  Same  as  in  Fig.  2  but  for  Category  4  lopen-cel led 
cumulus ) . 
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PRESSURE  (MB) 
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FIG.  7. 
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Same  as  in  Fig.  2  but  for  Category  6  (altocumulus 


PRESSURE  (MB) 


PRESSURE  (MB) 


PRESSURE  (MB) 


FIG.  10.  Same  as  in  Fig.  2  but  for  Category  9  (thick  cirrus). 
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1 1 .  Same  as  in  Fig . 
and  high  clouds) . 


2  but  for  Category  10  (chaotic  low, 


PRESSURE  (MB) 


CATEGORY  1  1 


SIGNIFICANCE 


FIG.  12.  Same  as  in  Fig.  2  but  for  Category  11  (organized 
low,  mid.  and  high  clouds). 


CATEGORY  12 


CATEGORY  1 2 


on  the  figures)  means  that  the  probability  is  5  percent  that  the 
given  rank-sum  discrepancy  between  the  category  and  the  ensemble 
is  due  to  chance.  The  smaller  the  significance  value  the  more 
significant  the  difference  between  the  category  and  the  ensemble. 

In  all  but  three  of  the  categories  (categories  6,  8,  and  12) 
at  least  half  of  the  levels  showed  a  significance  of  0.05  or 
better.  Category  6  (disorganized  mid  clouds)  was  the  most  infre¬ 
quent  category  with  only  seven  soundings  used  in  the  composite. 
Due  to  this  small  sample  the  composite  and  the  significance 
values  may  not  be  representative.  (The  same  may  also  be  true  of 
category  7  [organized  mid  clouds]  which  had  only  eleven  sound¬ 
ings).  Although  category  8  (thin,  high  clouds)  has  little 
significance  from  1000-550  mb,  it  does  have  considerable  signifi¬ 
cance  from  500-350  mb.  This  distribution  makes  sense  physically 
since  the  high  clouds  would  be  found  somewhere  within  these 
upper  levels.  It  also  makes  sense  that  category  12  (undecided) 
has  very  little  significance  since  it  is  a  composite  of  all  those 
cases  in  which  the  analyst  could  not  determine  a  category. 

The  significance  values  in  Figs.  2-13  show  that  the  compos¬ 
ites  have  significant  differences  when  compared  with  the  rest  of 
the  sample.  Significance  values  were  also  computed  to  determine 
the  significance  of  one  composite  when  compared  to  another  com¬ 
posite.  Due  to  the  large  number  of  comparisons  (65)  the  values 
are  not  presented  graphically  level  by  level.  Instead  values  are 


presented  in  Table  4.  Two  values  are  given  for  each  comparison. 
The  first  value  is  the  significance  value  for  the  comparison  as  a 
whole.  This  value  was  computed  by  averaging  the  test  scores  for 
each  of  the  15  pressure  levels  and  determining  the  significance 
for  this  average  score.  The  second  value  (in  parenthesis)  is  the 
number  of  pressure  levels  with  a  significance  of  0.05  or  better. 
For  example,  the  comparison  of  category  1  (clear)  versus  category 
2  (stratus)  shows  a  signicance  value  of  0.1235  for  the  comparison 
as  a  whole  (the  average  of  all  15  levels),  and  that  4  of  the  15 
levels  have  significance  values  of  0.05  or  better.  Only  13  of 
the  66  comparisons  had  fewer  than  five  levels  with  significance 
values  of  0.05  or  better. 

In  general  most  of  the  comparisons  show  a  considerable 
amount  of  significance  when  compared  to  each  other.  Two  of  the 
comparisons  which  showed  very  little  significant  difference  were 
2  (stratus)  versus  3  (stratocumulus) ,  and  9  (thick  cirrus)  versus 
11  (thick,  organized  multi-level  clouds).  The  similarity  in  2 
and  3  is  understandable  since  they  are  so  closely  related  physi¬ 
cally.  The  similarity  between  9  and  11  is  not  so  easily  under¬ 
stood.  It  is  most  likely  the  result  of  errors  in  the  satellite 
image  analysis,  i.e.  analyzing  multi-layered  cases  as  thick  cir¬ 
rus.  This  possibility  illustrates  the  limitations  of  both  the 
bogusing  technique  and  its  development  when  based  on  subjective 
satellite  image  analysis. 
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Table  4. 

Signi l icance 

(Number  of  Levels  with  0.05  or  Better) 


1  ” 

— 

\ 

1  ' 

2  > 

3  ! 

4  ! 

5  » 

6  ! 

7  ! 

8  ! 

9  !  10  i  11  ' 

1 

12351 

I 

I 

I 

» 

1 

I 

1 

!  l  i 

2 

1 

(4)  ' 

I 

1 

1 

i 

I 

I 

I 

!  1  ' 

1 

1007  i 

1912  ' 

1 

I 

i 

I 

I 

I 

3 

1 

(6)  ! 

(2)  ' 

! 

! 

i 

1 

l 

l 

!  !  i 

1 

0763  ' 

0173  ' 

.0219! 

I 

i 

l 

1 

1 

!  I  i 

4 

1 

(5)  ' 

(6)  ' 

(8)  ! 

l 

i 

i 

l 

l 

1  i  i 

1 

0068  ' 

0006  I 
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4.  DISCUSSION  AND  SUGGESTIONS 


\ 


Although  the  composites  are  statistically  significant,  the 
question  remains  as  to  whether  they  are  physically  significant, 
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and  if  they  would  produce  significantly  different  results  if 
implemented  in  a  numerical  model  initial  analysis.  The  second 
question  (model  impact)  will  not  be  addressed  here,  but  should  be 
explored  in  the  future.  The  first  question  (physical  signifi¬ 
cance)  will  be  discussed.  The  trend  from  dry  to  moist  and  the 
vertical  levels  with  relative  humidity  maxima  seem  to  qualita¬ 
tively  fit  the  worded  descriptions.  But  the  moistest  relative 
humidity  value  in  any  of  the  composites  is  92  percent  and  that 
occurs  only  once.  The  next  three  highest  values  are  only  90 
percent.  Obviously  in  the  real  atmosphere,  when  there  are 
significant  cloud  layers  present,  relative  humidity  at  or  near 
100  percent  would  be  observed. 

Four  possible  explanations  for  the  lack  of  relative  humidity 
values  near  saturation  in  the  composites  are:  1)  error  in  the 

satellite  image  analysis,  2)  inherent  spatial  variability  of 
moisture  in  the  atmosphere,  3)  radiosonde  measurement  error,  and 
4)  underestimated  relative  humidities  because  saturation  vapor 
pressure  for  ice  is  less  than  it  is  over  liquid  water.  Errors 
due  to  3)  and  4)  may  be  present  in  any  radiosonde  measurement 
and  thus  even  in  model  initial  analyses  derived  in  a  data  rich 
region.  Their  inclusion  in  the  bogusing  composites  does  not  add 
any  additional  error  to  the  analyses  when  compared  to  areas  using 
real-time  radiosonde  data. 

In  regards  to  the  first  explanation,  satellite  image  analy- 
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sis  error,  there  are  many  factors  which  could  contribute  to  this 
type  of  error.  The  first  is  simply  the  subjective  nature  of  the 
analysis.  For  example,  in  Garand  (1988)  ,  two  satellite  image 
analysts,  considered  to  be  experts  in  the  field,  agreed  only  37 
percent  of  the  time  on  a  strict  basis,  and  only  55  percent  of  the 
time  if  their  second  choice  was  considered  when  using  a  20-cate¬ 
gory  cloud  classification  scheme  (a  study  is  currently  underway 
at  Point  Mugu  to  study  the  differences  between  analysts  using  the 
categorization  scheme  presented  in  this  paper) . 

Another  contributing  factor  in  this  type  of  error  is  the 
quality  of  the  satellite  imagery  used.  As  stated  earlier,  NOAA-2 
polar  orbiting  satellite  infrared  and  visual  image  mosaics  for 
the  Northern  Hemisphere  were  used  in  this  study.  Unfortunately, 
the  resolution  of  the  imagery  was  somewhat  limiting  due  to  the 
large  area  contained  in  the  image.  Also,  the  satellite  pass  over 
the  sounding  point  was  as  much  as  six  hours  off  from  the  sounding 
time.  In  general  only  those  soundings  taken  within  three  hours 
of  the  satellite  pass  were  used.  The  exception  was  when  very 
little  synoptic  change  was  observed  between  images.  In  some 
cases  even  when  the  soundings  were  within  two  to  three  hours  of 
the  satellite  pass,  speculation  had  to  be  made  on  the  movement, 
development,  or  dissipation  of  cloud  features  near  the  sounding 
point.  This  was  particularly  true  in  the  case  of  narrow  frontal 
zones  in  the  vicinity  of  the  sounding  point.  A  new  data  base  of 
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higher  quality  imagery  taken  within  one  hour  of  the  sounding  time 
could  significantly  reduce  this  type  of  error.  Efforts  are 
currently  underway  to  obtain  such  a  data  set. 

The  second  explanation,  inherent  spatial  variablity,  is 
almost  certainly  a  reason  for  the  lack  of  relative  humidities 
near  saturation  in  the  composites.  This  variability  manifests 
itself  in  two  ways.  First,  cloud  patterns  do  not  always  fit  into 
a  finite  number  of  categories.  There  may  be  enough  similari  ty 
from  one  case  to  the  next  to  warrant  having  distinct  groupings, 
but  there  will  always  be  cases  that  do  not  fit  the  categoriza¬ 
tion.  For  example,  56  of  the  471  soundings  used  in  this  study 
were  considered  unidentifiable  under  the  scheme  presented  in  this 
paper.  Lyons  (1986b),  in  determining  composites  for  the  North¬ 
east  Pacific,  tried  to  overcome  this  problem  by  selecting  only 
'classical'  examples  for  each  of  his  categories.  As  a  result, 
his  composites  were  based  on  only  two  to  ten  soundings,  thus 
leaving  open  the  question  of  their  validity  in  representing  the 
varied  conditions  found  in  the  atmosphere. 

The  second  manifestation  of  spatial  variability  is  in  the 
vertical.  Even  within  easily  identifiable  categories  there  will 
be  vertical  variations  in  the  height  and  thickness  of  the  cloud 
layers.  This  variation  produces  thicker  but  drier  moist  layers 
in  the  composites.  The  total  integrated  moisture  within  the 
column  of  air  may  be  similar  for  the  composite  and  the  individual 
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cases,  but  no  cloud  layer  would  be  identifiable  in  the  composite. 
This  could  present  a  problem  in  a  numerical  model  where  parame- 
terizations  of  precipitation  and  radiation  processes  are  not  used 
until  the  relative  humidity  nears  saturation.  One  possible  sug¬ 
gestion  to  overcome  this  problem  of  vertical  smoothing  would  be 
to  determine  the  distribution  of  cloud  height  and  thickness  based 
on  the  sounding  relative  humidities  and  create  a  saturated  cloud 
layer  at  the  most  frequently  occurring  levels  with  an  average 
thickness  taken  from  the  data.  Moisture  would  then  have  to  be 
taken  out  in  other  areas  in  order  to  preserve  the  same  amount  of 
total  moisture  content  within  the  column.  This  technique  will  be 
investigated  more  fully  in  the  future. 

5 .  SUMMARY 

Vertical  moisture  profiles  for  use  in  satellite  bogusing  of 
moisture  into  numerical  model  initial  analyses  were  presented  for 
the  North  Atlantic  Ocean  region.  The  profiles  were  based  on  a  12 
category  cloud  classification  scheme  and  computed  from  471  sound¬ 
ings  taken  in  the  North  Atlantic  during  January,  July,  and 
October  of  1974.  Assignment  of  soundings  to  one  of  the  cloud 
categories  was  done  by  subjective  analysis  of  Northern  Hemisphere 
infrared  and  visual  satellite  image  mosaics  taken  from  the  NOAA-2 
polar  orbiting  satellite.  Composites  for  the  categories  and 
results  of  statistical  significance  tests  were  presented.  The 
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composites  appeared  to  vary  appropriately  according  to  their  dif¬ 
ferent  descriptions.  The  significance  tests  showed  that  the  vast 
majority  of  the  composites  were  statistically  significant  when 
compared  to  the  sample  and  when  compared  to  each  other. 

A  lack  of  relative  humidity  values  near  saturation  was 
considered  cause  for  concern.  Errors  in  satellite  image  analy¬ 
sis,  inherent  spatial  variability  of  moisture  in  the  atmosphere, 
radiosonde  measurement  error,  and  differences  in  saturation  vapor 
pressure  over  liquid  water  and  ice  were  considered  as  possible 
reasons  for  the  lack  of  relative  humidity  values  near  saturation 
in  the  composites.  More  highly  resolved  satellite  imagery  closer 
to  the  sounding  times  was  suggested  as  a  possible  way  to  signifi¬ 
cantly  decrease  this  error.  A  technique  for  retaining  cloud 
layers  in  the  compositing  process  was  suggested  as  a  possible 
means  of  reducing  the  smoothing  of  relative  humidity  in  the 
vertical  due  to  vertical  moisture  variability. 

ACKNOWLEDGEMENTS 

The  author  would  like  to  thank  all  those  who  assisted  in 
this  work:  Paul  Tag  and  Allan  Caughey  of  NEPRF  for  supplying 
much  of  the  sounding  data;  Jay  Rosenthal  for  developing  the  new 
cloud  classification  scheme;  Chuck  Fisk  for  providing  computer 
code  and  expertise  on  nonparametr ic  statistics;  and  Roger  Helvey, 
Alan  Fox,  and  Richard  Szymber  for  their  discussions  and  ideas. 


27 


REFERENCES 


Chu,  R. ,  1977:  Humidity  analyses  for  operational  prediction 

models  at  the  National  Meteorological  Center.  NMC  Office 
Note  140,  Part  I,  14  pp . 

Qarand ,  L.,  1988:  Automated  recognition  of  oceanic  cloud 

patterns.  Part  I:  Methodology  and  application  to  cloud 
climatology.  J.  Climate,  1,  20-39. 

Lyons,  S.  W. ,  1986a:  Satellite  derived  moisture  fields  as  input 
to  operational  numerical  forecast  models:  NMC  case  study 
April  5-6,  1986.  Geophysical  Sciences  Tech.  Note  No.  118, 
Pacific  Missile  Test  Center,  Point  Mugu,  CA  93042-5000,  14 

pp. 

- ,  1986b:  Satellite  derived  moisture  fields  as  input  to 

operational  numerical  forecast  models:  A  modification  of 
NMC  relative  humidity  profiles.  Geophysical  Sciences  Tech. 
Note  No.  122,  Pacific  Missile  Test  Center,  Point  Mugu,  CA 
93042-5000,  6  pp . 

Smigielski,  F.,  T.  Burtt  and  S.  Hirsch,  1982:  Moisture  bogus 
program.  N0AA/NESS/SAB  Office  Operations  Note  (Chapter 
225) ,  22  pp. 

Timchalk,  Andre,  1986:  Satellite-derived  moisture  profiles. 

N0AA  Tech.  Report  NESDIS  24,  U.S.  Dept,  of  Commerce,  N0AA, 
Rockville,  MD ,  20852,  60  pp . 


DISTRIBUTION 


ASST.  FOR  ENV.  SCIENCES 
ASST.  SEC.  OF  THE  NAVY  (RAD) 
ROOM  5E731,  THE  PENTAGON 
WASHINGTON,  DC  20350 


OFFICE  OF  NAVAL  RESEARCH 
ENV.  SCI.  PROGRAM,  CODE  112 
ARLINGTON,  VA  22217-5000 


COMMANDING  OFFICER 
NAVAL  RESEARCH  LAE 
ATTN:  LIBRARY,  CODE  2620 
WASHINGTON,  DC  20390 


COMNAVOCEANCOM 
ATTN:  CODE  N5 
JCSSC,  MS  39529-5000 


SUPERINTENDENT 
LIBRARY  REPORTS 
U.S.  NAVAL  ACADEMY 
ANNAPOLIS,  MD  21402 


LIBRARY 

NAVAL  POSTGRADUATE  SCHOOL 
MONTEREY,  CA  93943-5002 


COMMANDER,  SPAWARSYSCOM 
ATTN:  CODE  141 
NAT.  CTR.  #1 

WASHINGTON,  DC  20363-5100 


AFGWC/DAPL 

OFFUTT  AFB,  NE  68113 


CHIEF  OF  NAVAL  RESEARCH 
LIBRARY,  CODE  012321 
BALLSTON  TONER  #1 
800  QUINCY  ST. 

ARLINGTON,  VA  22217-5000 


OFFICE  OF  NAVAL  RESEARCH 
ATTN:  HEAD,  OCEAN  SCIENCES  DIV 
CODE  1122 

ARLINGTON,  VA  22217-5000 


COMMANDING  OFFICER 
NORDA,  CODE  335 
JCSSC,  MS  39529-5004 


COMMANDING  OFFICER 
NAVAL  OCEANOGRAPHIC  OFFICE 
JCSSC,  MS  39522-5001 


CHAIRMAN 

OCEANOGRAPHY  DEPT. 
U.S.  NAVAL  ACADEMY 
ANNAPOLIS,  MD  21402 


COMMANDER 

NAVAIRSYSCCM,  CODE  526W 
WASHINGTON,  DC  20361-0001 


COMMANDER 
PACMISTESTCEN 
GEOPHYSICS  OFFICER 
PT.  MUGU,  CA  93042 


AFGL/LY 

HANSCOM  AFB,  MA  01731 


OFFICE  OF  NAVAL  RESEARCH 
CODE  1122AT,  ATMOS.  SCIENCES 
ARLINGTON,  VA  22217-5000 


OFFICE  OF  NAVAL  TECHNOLOGY 
ONR  (CODE  22) 

800  N.  QUINCY  ST. 
ARLINGTON,  VA  22217-5000 


COMMANDER 

NAVAL  OCEANOGRAPHY  COMMAND 
JCSSC,  MS  39529 


COMMANDING  OFFICER 
FLENUMOCEANCEN 
MONTEREY,  CA  93943-5005 


NAVAL  POSTGRADUATE  SCHOOL 
METEOROLOGY  DEPT. 
MONTEREY,  CA  93943-5000 


COMMANDER,  SPAWARSYSCOM 
ATTN:  CODE  PMW  141-Bl 
WASHINGTON,  DC  20363-5100 


USAFETAC/TS 
SCOTT  AFB,  IL  62225 


COMMANDER/DIRECTOR 
ATMOS.  SCI.  LABORATORY 
ATTN:  SLCAS-AE-A  (HANSEN) 
WHITE  SANDS  MR,  NM  88002-5501 


DIRECTOR  (6) 

DEFENSE  TECH.  INFORMATION 
CENTER,  CAMERON  STATION 
ALEXANDRIA,  VA  22314 


DIRECTOR 

NATIONAL  EARTH  SAT.  SERV/SEL 
FB-4,  S321B 
SUITLAND,  MD  20233 


DIRECTOR 

ATLANTIC  MARINE  CENTER 
COAST  &  GEODETIC  SURVEY,  NOAA 
439  W.  YORK  ST. 

NORFOLK,  VA  23510 


DIRECTOR 

TECHNIQUES  DEVELOPMENT  LAB 
GRAMAX  BLDG. 

8060  13TH  ST. 

SILVER  SPRING,  MD  20910 


SCRIPPS  INSTITUTION  OF 
OCEANOGRAPHY,  LIBRARY 
DOCUMENTS/REPORTS  SECTION 
LA  JOLLA,  CA  92037 


CHAIRMAN,  METEOROLOGY  DEPT. 
PENNSYLVANIA  STATE  UNIV. 

503  DEIKE  BLDG. 

UNIVERSITY  PARK,  PA  16802 


LIBRARY,  AUSTRALIAN  NUMERICAL 
METEOROLOGY  RESEARCH  CENTER 
P.O.  BOX  5089A 
MELBOURNE,  VICTORIA,  3001 
AUSTRALIA 


NOAA-NESDIS  LIAISON 
ATTN:  CODE  SC2 
NASA -JOHNSON  SPACE  CENTER 
HOUSTON,  TX  77058 


CHIEF 

MARINE  &  EARTH  SCI.  LIBRARY 
NOAA,  DEPT.  OF  COMMERCE 
ROCKVILLE,  MD  20852 


DIRECTOR 

GEOPHYS.  FLUID  DYNAMICS  LAB 
NOAA,  PRINCETON  UNIVERSITY 
P.O.  BOX  308 
PRINCETON,  NJ  08540 


LIBRARY  ACQUISITIONS 
NCAR,  P.O  BOX  3000 
BOULDER,  CO  80307 


COLORADO  STATE  UNIVERSITY 
ATMOSPHERIC  SCIENCES  DEPT. 
ATTN:  LIBRARIAN 
FT.  COLLINS,  CO  80523 


DR.  CLIFFORD  MASS 
DEPT.  OF  ATMOS.  SCIENCES 
UNIV.  OF  WASHINGTON 
SEATTLE,  WA  98195 


EUROPEAN  CENTRE  FOR  MEDIUM 
RANGE  WEATHER  FORECASTS 
SHINFIELD  PARK,  READING 
BERKSHIRE  RG29AX,  ENGLAND 


DIRECTOR 

NATIONAL  METEORO.  CENTER 
NWS,  NOAA 
WWB  W32,  RM  204 
WASHINGTON,  DC  20233 


NATIONAL  WEATHER  SERVICE 
WORLD  WEATHER  BLDG.,  RM  307 
5200  AUTH  ROAD 
CAMP  SPRINGS,  MD  20023 


CHIEF 

MESOSCALE  APPLICATIONS  BRANCH 
NATIONAL  EARTH  SAT.  SERV. 

1225  W.  DAYTON 
MADISON,  WI  53562 


HEAD,  ATMOS.  SCIENCES  DIV. 
NATIONAL  SCIENCE  FOUNDATION 
1800  G  STREET,  NW 
WASHINGTON,  DC  20550 


UNIVERSITY  OF  WASHINGTON 
ATMOSPHERIC  SCIENCES  DEPT. 
SEATTLE,  WA  98195 


AMERICAN  METEORO.  SOCIETY 
METEOR.  &  GEOASTRO.  ABSTRACTS 
P.O.  BOX  1736 
WASHINGTON,  DC  20013 


30 


